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Abstract 

The Pierre Auger Observatory is a detector for ultra- high energy cosmic rays. It 
consists of a surface array to measure secondary particles at ground level and a fluo- 
rescence detector to measure the development of air showers in the atmosphere above 
the array. The "hybrid" detection mode combines the information from the two sub- 
systems. We describe the determination of the hybrid exposure for events observed 
by the fluorescence telescopes in coincidence with at least one water-Cherenkov de- 
tector of the surface array. A detailed knowledge of the time dependence of the 
detection operations is crucial for an accurate evaluation of the exposure. We dis- 
cuss the relevance of monitoring data collected during operations, such as the status 
of the fluorescence detector, background light and atmospheric conditions, that are 
used in both simulation and reconstruction. 

Key words: Ultra high energy cosmic rays, Pierre Auger Observatory, Extensive 
air showers, Trigger, Exposure, Fluorescence detector, Hybrid 



1 Introduction 



The Pierre Auger Observatory has been designed to investigate the origin and 
the nature of Ultra High Energy Cosmic Rays. It consists of a large array 
of about 1600 surface stations (the SD array) covering an area of 3000 km 2 
for detecting the secondary particles of the air shower at ground level by 
means of the Cherenkov radiation they produce in water. The ground array 
is overlooked by 24 air fluorescence telescopes (the FD system), grouped in 
4 enclosures each consisting of 6 optical telescopes. These devices are used 
to observe the longitudinal profile of cosmic ray showers on clear moonless 
nights. The Observatory, located outside the town of Malargue, in the Province 
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of Mendoza, Argentina, has been taking data stably since January 2004 while 
the construction was proceeding. The construction was completed in mid 2008. 
Details of the design, construction and performance of the Observatory can 
be found in [1,2,3]. 

The Auger detector has been conceived with a cross-triggering capability. Data 
are retrieved from both detectors whenever either system is triggered 1 1. The 
surface array and the fluorescence telescopes allow the reconstruction of ex- 
tensive air showers with two independent measurements. The combination of 
information from the two detection subsystems enhances the reconstruction 
capability with respect to the individual detector components [5,6]. This tech- 
nique is called "hybrid" detection and the determination of the exposure of 
the Observatory under this mode is the subject of the present paper. The data 
period used for this purpose is between November 2005 and May 2008. The ex- 
posure calculated here is the same used for the energy spectrum measurement 
published in ([7]). 

The paper is organized as follows. In section 2 we describe the hybrid detection 
method. Section 3 addresses the energy spectrum and the relevance of the 
hybrid exposure to its determination. The effective data taking time in the 
hybrid detection mode, i.e the hybrid on-time, and the different components 
contributing to it are discussed in section 4. The Monte Carlo simulation used 
for the evaluation of the hybrid exposure is described in section 5. In section 
6 we describe the event selection and make comparisons with data to validate 
the Monte Carlo simulation. Finally in section 7 we show the hybrid exposure 
as a function of primary energy and in section 8 we summarize. 



2 Hybrid data analysis 

A hybrid event is an air shower that is simultaneously detected by the fluores- 
cence detector and the surface array. If an air shower independently triggers 
both detectors the event is tagged as a golden hybrid and these events can be 
fully reconstructed in both detection modes. In the SD the energy density of 
shower particles at ground level is used to determine the cosmic ray energy. 
In the FD the observation of the longitudinal profile of the shower allows the 
measurement of the calorimetric energy of the primary particle. This event 
sample, though small with respect to the SD sample, is very important since 
it constitutes the base data set for the energy calibration of the SD events 
[8,6]. 



1 Details about the triggers implemented can be found in [4] for the SD and in [3] 
for the FD. 
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The fluorescence detector, having a lower energy threshold, may promote a 
sub-threshold trigger in the SD. In this case, surface stations are matched 
by timing and location even though they do not fulfil the conditions for an 
independent SD trigger. This is an important capability because these sub- 
threshold hybrid events would not have triggered the array otherwise. Here 
the energy reconstruction relies uniquely on the calorimetric energy from the 
longitudinal profile. 

Like golden hybrids these events suffer statistical limitations, but they are of 
particular interest because they allow an extension of the measurement of the 
energy spectrum into a region where the SD is not fully efficient [4]. They have 
superior qualities with respect to "monocular" FD events (those without SD 
information), because of the precise measurement of the shower geometry [3]. 

In the FD, cosmic ray showers are detected as a sequence of triggered pixels 
in a matrix of photomultipliers. This sequence allows the determination of 
the shower-detector plane (SDP), the plane that includes the location of the 
fluorescence detector and the line of the shower axis, with a typical uncertainty 
of the order of a few tenths of a degree. Then the determination of the shower 
geometry relies on the arrival times of photons in the individual pixels [3]. 
In the monocular reconstruction the accuracy degrades when the measured 
angular speed does not change significantly over the observed track length. In 
such cases the shower axis can be largely under-determined within the SDP, 
thus giving large uncertainties in the reconstruction of the arrival direction 
and the impact point at ground level. This further leads to uncertainties in 
other shower parameters and in particular in the reconstructed shower energy. 

The hybrid approach supplements the traditional FD direction fitting method 
with the arrival time of the shower at the ground measured by a single SD 
station. This results in a remarkable improvement in the determination of the 
shower geometry, as illustrated in figure 1 where the impact points at ground 
level corresponding to mono and hybrid reconstruction methods are shown for 
a typical event. Accurate knowledge of the shower arrival time at ground level 
removes a degeneracy in the traditional FD monocular approach that uses 
pixel timing to reconstruct the shower axis. In hybrid mode, the resolution of 
the direction and of the position of the impact point at the ground are better 
than 0.6° and 50 m respectively [9,10,5]. 

The total energy of each event is obtained by combining the knowledge of 
the detector response with monitoring data describing the atmospheric con- 
ditions [3]. Once the geometry is known, the observed energy deposit profile 
is reconstructed taking into account the scattering and the absorption of light 
during its propagation in the atmosphere and the presence of forward-emitted 
and scattered Cherenkov light. The method used is described in detail in [11]. 
The energy released in the electromagnetic part of the air shower is estimated 



S 




3 4 5 6 7 8 

x [km] 

Fig. 1. Determination of the impact point at the ground for a single event using both 
the mono and hybrid reconstruction methods. The event has been detected by the 
Los Morados FD site: the downward-going arrow points towards the direction of the 
site and the two lines show the uncertainty of the SDP plane at ground level. The 
small (long) elongated ellipse represents the uncertainty on the core position in the 
hybrid (mono) reconstruction. The arrows indicate the reconstructed directions in 
the two cases, their length being proportional to the sine of the reconstructed zenith 
angle. The open (full) circles show the active (triggered) SD stations. Triggered 
stations are shown with a radius proportional to the logarithm of the signal. 

by fitting a Gaisser-Hillas function [12] to the reconstructed energy deposit 
profile and integrating it over the entire range of atmospheric depth. Finally, 
the total energy of a shower is derived after correcting for the invisible energy 
carried away by neutrinos and high energy muons [13]. After quality selection, 
the energy resolution (defined as an event-to-event statistical uncertainty) of 
the fluorescence detector is better than 10% [5]. This value has been calcu- 
lated from simulations and has been cross checked with hybrid-stereo events, 
i.e. those events which are detected and reconstructed in hybrid mode by more 
than one FD-site[3]. The energy resolution turned out to be energy indepen- 
dent in the whole range. 
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Table 1 

Current estimates of the systematic uncertainties affecting energy reconstruction. 
Values from [6]. 





uncertainty % 




uncertainty % 


fluorescence yield (FY) 


14 


quenching effect on FY 


5 


FD absolute calibration 


9 


FD wavelength response 


3 


molecular attenuation 


1 


aerosol attenuation 


7 


multiple scattering model 


1 


FD reconstruction method 


10 


invisible energy 


4 






total 






22 



Systematic uncertainties in the energy determination are related to the detec- 
tor, to the atmosphere and to the reconstruction procedure. They are sum- 
marized in Tab. 1. All these uncertainties are found to be independent. A 
total uncertainty of about 22% [6] is estimated by summing the individual 
contributions in quadrature. 



3 Energy Spectrum with hybrid events 



The aperture of a cosmic ray instrument is per se a figure of merit of its 
observation capability. The time integrated aperture is commonly referred 
to as the exposure. In this section we discuss the relevance of the exposure 
for the energy spectrum measurement. This is of particular concern in the 
case of a detection based on fluorescence, such as the hybrid case, where the 
time variations of the detection and the inherent energy dependence make an 
accurate determination of the exposure a key task. 

The flux of cosmic rays J as a function of energy is defined as: 

J(E) = d ^ AN scl (E) 1 . 
v ; dEdAdttdt AE £(EY V ; 

where iVi nc is the number of cosmic rays with energy between E and E+dE in- 
cident on a surface element dA, within a solid angle dQ and time dt. AN SC \(E) 
is the number of detected events passing the selection criteria in the energy bin 
centered around E, having width AE. S(E) represents the energy- dependent 
exposure of the detector at the same selection level. 

The exposure, as a function of the energy of primary particle, can be written 
as: 

S(E) = I I I e(E t t 7 9,<f>,x,y) cos6 dS dfi dt = [ A(E,t) dt; (2) 

JT Jn JSgen JT 
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where e is the detection efficiency including the different steps of the analysis, 
i.e trigger, reconstruction and quality cuts, and dS = dxxdy is the horizontal 
surface element, dfl = sin#d#d0 and Q are respectively the differential and 
total solid angles. The generation area S gen has been chosen large enough to 
exclude any possible event detection and reconstruction outside it. A(E,t) 
is the instantaneous aperture of the detector which depends on the detector 
configuration at the time t. 

The detector configurations of the Observatory have been continuously chang- 
ing over the period of data collection for the hybrid spectrum. As construction 
of the SD progressed, the number of stations in operation increased. Further- 
more, even in a steady configuration, some SD stations are temporarily out of 
service at any one time. The SD status is monitored by updating each second 
the list of "active" stations. In principle the change in SD configuration is 
straightforward to handle since the aperture is proportional to a geometric 
area. In the case of a single missing SD station, the effective area is slightly 
changed by about 2 km 2 at full efficiency [4] . 

The FD detector configuration also changed with time during the construction 
phase, with the number of telescopes changing from 12 to 24. In addition, a 
correction ring lens was added to each telescope during the first two years of 
data taking. Thus, parts of the data have been collected with different optical 
configurations. During nightly operations individual telescopes are sometimes 
deactivated because of increasing sky brightness, bad weather conditions or 
hardware failures. Finally, the FD response is influenced by atmospheric con- 
ditions such as the concentration of aerosols and cloud coverage. 

To properly take into account all the detector configurations and their time 
variability a sample of events which reproduce the exact conditions of the 
experiment (i.e its actual sequence of configurations and on-time) has been 
simulated. This method, referred to as Time Dependent Detector Simulation, is 
described in the next sections. Given a set of N simulated events generated on 
an area S gen within the time interval T, the exposure eq. (2) can be calculated 
numerically via 

S(E iec ) = 2tt S gcn T V ^^ rec,CQS y. cos 6); Acos^; (3) 



where n denotes the number of events that fulfill the selection criteria de- 
scribed in Sec. 6. The exposure is calculated as a function of reconstructed 
energy, E Tec , to correct for distortions of the steep energy spectrum due to the 
finite resolution of the energy reconstruction (see e.g. [14,15] and Sec. 6.3). 
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4 Hybrid on-time 



The efficiency of fluorescence and hybrid data taking is influenced by many 
effects. These can be external, e.g. lightning or storms, or internal to the data 
taking itself, e.g. DAQ failures. For the determination of the on-time of the 
Pierre Auger Observatory in the hybrid detection mode it is therefore crucial 
to take into account all these occurrences and derive a solid description of the 
data taking time sequence. 

Data losses and inefficiencies can occur on different levels, from the smallest 
unit of the FD, i.e. one single photomultiplier (pixel) readout channel, up to 
the highest level, i.e. the combined SD-FD data taking of the Observatory. To 
perform the time dependent detector simulation we have to take into account 
all known disturbances and then derive the on-time of the hybrid detection 
mode. To achieve this aim we rely on a variety of monitoring information 
and the data set itself. As a compromise between accuracy and stability we 
derived the complete detector status down to the single pixel for time intervals 
Tbm = 10 min. 



4-1 Telescope dependent sources 

The active time of FD data acquisition is calculated using a minimum bias 
data stream with a less restrictive trigger condition. This data file includes 
sub-threshold FD events and is recorded at an event rate about 8 times higher 
than the standard rate of about 1 event per FD-site per minute. 

Even if the DAQ is running, the shutters of the telescope might be closed 
due to bad weather alarms from the slow control system or other failsafe 
mechanisms. To determine the status of the shutters we use the information on 
night sky background level provided by algorithms implemented in the front- 
end electronics boards. Every 30 s data from each PMT channel is written 
to a monitoring data file which records parameters including ADC-variance, 
baseline, First Level Trigger (FLT) threshold and trigger frequency for each 
pixel [3]. 

The ADC- variance distribution from these data is shown in figure 2. Back- 
ground data is also collected during the nightly relative calibration runs, i.e. 
with closed shutters (see the upper- right panel in figure 2). A mean value of 
about 3.5 ADC 2 is obtained in these condition^. For each time interval the 
efficiency of open shutters is then derived as: 



Muons hitting the pixel camera is the main source of the noise triggers. 
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Fig. 2. Distribution of background variances. The main contribution to this back- 
ground noise is the night-sky background light coming from stars and the direct and 
scattered moonlight. The upper-right panel shows a magnified view of the low vari- 
ance region superimposed on data recorded with closed shutters (shaded histogram). 
The arrow shows the variance threshold used to select good data. 
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where T open denotes the time (for a given telescope) for which the mean vari- 
ance over the whole camera is larger than 8 ADC 2 . If background data are 
not available, no efficiency is calculated. The status flag 5 tc i is then set to 0. 

The deadtime due to the finite readout speed of the DAQ system must also 
be taken into account. The deadtime is stored on an event-by-event basis in 
the output of the FD data acquisition. For each telescope, this deadtime 
is converted into an efficiency of detecting cosmic ray data in the given time 
interval by: 



T^dead 

^DAQ = 1 - ^ (5) 
J- DAQ 



where Tdaq is the total running time of the DAQ in the given time interval. 
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4-2 FD-site dependent sources 



Currently two possible sources of inefficiency are known to affect the data 
taking at the FD-site level. 

The first is due to the atmospheric monitoring system. An FD veto is set by the 
Lidar system before performing laser shots in the field of view of a fluorescence 
detector. The cumulative Lidar veto time is stored on an event-by-event basis 
in the data files. This deadtime T^^ r is converted into an efficiency by: 



ydead 

^Lidar = 1 - TfT^- (6) 
-*DAQ 

This efficiency can be interpreted as the probability of a cosmic ray event 
falling outside the Lidar vetoed period. 

To extend the hybrid detection capability below the SD trigger threshold [4,16], 
all the FD triggers are sent to and processed by the central data acquisition 
system (CDAS). It reads out the portion of the surface array closest to the 
relevant fluorescence building. Then FD and SD data streams are merged to 
form hybrid events. A source of inefficiency comes from the protection algo- 
rithm implemented in the CDAS to prevent the acquisition of long periods of 
excessive event ratedjll This veto mechanism induces the loss of hybrid events. 
An estimate of the event loss probability in a given time interval is calculated 
by comparing events from the FD data files and from the final merged hybrid 
files (which only include those sent to CDAS). This recovery mechanism is 
energy dependent as it is related to the SD trigger probability [17] and is ac- 
counted for on an average basis. (eT 3vcto (s, t)) is the resulting average efficiency 
for each FD site s and time t. 



4.3 CDAS status 



CDAS inefficiencies must also be taken into account. The surface detector 
array is constantly monitored and a very detailed description of the array 
status is available with a time resolution of 1 second. In addition to the usually 
very localized problems of single SD stations, time periods with trigger related 
problems [4] are excluded in the hybrid on-time via the CDAS status flag 

3 Lightning and other noise events may cause higher FD trigger rates which would 
cause significant deadtime for the surface array due to the finite readout time of the 
array. 
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^cdas- Given a constant rate of hybrid events A, the probability P that the 
time interval between two consecutive hybrid events is larger than T is given 
by P(T) = e~ XT . Taking A « 1.7 X 10" 2 Hz (1 event per minute) and T = 600 
sec, then P = 3.7 x 10~ 5 . Based on this calculation, an additional check is 
performed requiring at least one hybrid event per 10 min time interval. 

4-4 Results and cross-checks 

For each time t in a given time slot of duration Tbi n , the fraction of operational 
time f(i,t), for the telescope i belonging to the FD site s, can be written as: 

f(i, t) = £ sh uttcr(«, t) ■ e DA Q(i, t) ■ 5 te \(i, t) (7) 

'^Lidar 

(s,t) ■ (£T3vcto(s,£)) 
•^CDAs(^) 

where the e's identify the efficiencies due to the different sources and the <5's 
are status flags (5 = [0, 1]). All the expected sources of inefficiencies have been 
described in detail in the previous sections. 

The time evolution of the full hybrid duty-cycle over 3 years during the con- 
struction phase of the observatory is shown in figure 3. It shows the on-time 
fraction, defined as the ratio of the overall on-time to the time duration of each 
interval. To avoid pile-up effects in the plot, time bins are chosen to coincide 
with FD data-taking shifts. Data-taking is currently limited to dark periods 
with moon-fractions smaller than 60% as seen by each individual telescope: 
this leads to about 16 nights of data taking per moon-cycle. The scheduled 
data-taking time fraction is also shown in figure 3 (gray line). A seasonal 
modulation is clearly visible, since higher fractions are observed in the austral 
winter during which the nights are longer. 

Note that the FD-site at Los Morados became operational in May 2005 and 
that at Loma Amarilla started in March 2007. After the initial phase of com- 
missioning, the mean on-time is about 12% for all FD-sites, which corresponds 
roughly to about 70% of the scheduled time fraction. This efficiency is primar- 
ily due to weather effects with a minor part determined by detector effects. 

A validation of the on-time determination and an estimate of its systematic 
uncertainty has been performed using data from the Central Laser Facility 
(CLF) [18]. These data are embedded in the standard FD data stream. As 
CLF laser shots can be observed from all FD-sites, one can calculate the 
conditional probability of recording the laser signal in a particular site s given 
at least one other observation in any other site. The expected number of 
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Fig. 3. Time evolution of the average hybrid on-time fraction during the construc- 
tion phase of the Pierre Auger Observatory. Both the seasonal modulation and 
the starting of commissioning phases of the different FD-sites are visible. Gray 
line represents the scheduled data-taking time fraction limited to the nights with 
moon- fraction lower than 60%. 



laser shots in site s can be derived from the on-time of the telescope pointing 
to the CLF. The laser observation probability is obviously dependent on the 
transmission coefficient of the atmosphere. The probability of observing a laser 
during aerosol-free periods, i.e with vertical aerosol optical depth VAOD « 0, 
is expected to be 100%. A small deviation from this value of about 4% was 
found and the on-time has been corrected accordingly to account for possibly 
lost periods. 



5 Monte Carlo Simulation 



For the calculation of the hybrid exposure, the size of the simulated event sam- 
ple is crucial for acceptable statistical and systematic uncertainties. For this 
purpose the simulation activity followed a graded approach with full Monte 
Carlo analysis for specific studies, like the trigger efficiency, and fast simula- 
tions, validated with the full Monte Carlo method, when high statistics were 
required. 
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Fig. 4. Relative hybrid trigger efficiency from hybrid simulation for proton and iron 
primaries. The hybrid trigger efficiency calculated using data is also shown. 

5.1 Trigger efficiency 

A complete Monte Carlo hybrid simulation has been performed to study the 
trigger efficiency and the detector performance. The simulation sample con- 
sists of about 6000 proton and 3000 iron CORSIKA [19] showers with en- 
ergies ranging between 10 17 and 10 19 ' 5 eV. These energies are of particular 
interest for the trigger studies since they cover both SD and hybrid thresh- 
olds. The showers have been generated using respectively QGSJET-II [20,21] 
and FLUKA[22] as high and low energy hadronic interaction models. The FD 
simulation chain [23] reproduces in detail all the physical processes involved 
in the fluorescence technique. It includes the generation of fluorescence and 
Cherenkov photons in the atmosphere, their propagation through the atmo- 
sphere to the telescope aperture, the ray-tracing of photons in the Schmidt 
optics of the telescopes, and the simulation of the response of the electronics 
and of the multi-level trigger. The surface detector response is simulated using 
Geant4 [24] within the framework provided by the Auger Offline software [25] . 
For this particular purpose we assume the SD array is fully operational and 
deployed. 

In figure 4 it is shown the hybrid trigger efficiency, i.e. the probability of 
detecting a fluorescence event in coincidence with at least one triggered SD 
station, is flat and equal to 1 at energies greater than 10 18 eV, independent 
of primary mass. The difference between proton and iron primaries increases 
at lower energies but is negligible at energies as low as 10 17 ' 5 eV. Protons 
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are slightly more efficient than iron primaries at the lowest energies. This is 
mainly due to the larger fraction of proton events interacting deeper in the 
atmosphere. The hybrid trigger efficiency from fluorescence data is also shown 
in figure 4. Only events landing on an active part of the surface detector have 
been selected and minimal quality cuts have been applied in order to have a 
reliable reconstructed energy and to safely derive the trigger probability curve. 
Data and simulation consistently show that a fluorescence event is always 
hybrid for energies larger than 10 18 eV. 

In addition, the probability of a shower triggering a given SD station has 
been studied as a function of primary cosmic ray energy, mass, direction and 
distance to the shower axis, and a set of "Lateral Trigger Probability" (LTP) 
functions have been derived and parameterised [26]. For a vertical proton 
primary shower, each station is on average fully efficient within a distance of 
750, 1000, 1300, and 1600 m at energies of 10 17 5 eV, 10 18 eV, 10 18 5 eV and 
10 19 eV, respectively. Details on this study are discussed in [26]. 



5.2 Fast simulation 

To follow and reproduce the time dependence of the hybrid exposure, each 
detector configuration must be taken into account. This approach requires a 
large number of simulations. The method used to achieve this goal within a 
reasonable computational time relies on the simulation of longitudinal shower 
profiles generated with CONEX [27], a fast generator based on CORSIKA 
[19] shower code. After the simulation of the first few ultra-high energy in- 
teractions, CONEX switches to numerical solutions of the underlying cascade 
equations that describe the evolution of the different shower components. Al- 
though this method is extremely fast, the most important features provided 
by full Monte Carlo simulations, including shower to shower fluctuations, are 
very well reproduced [27,28]. 

The simulation of the FD response proceeds as in the full method discussed 
above. Since no ground level particles are generated by CONEX, the SD re- 
sponse cannot be directly simulated. In this case the SD trigger is reproduced 
using the LTP parameterisation functions. The actual status of the SD array 
is retrieved using the time of each simulated event. The event trigger proba- 
bility is then calculated as the convolution of all the LTPs of the working SD 
stations. This is particularly important for low energy and inclined events. 

The SD timing information needed in the hybrid reconstruction mode is pro- 
vided by a simplified simulation (i.e. SdSimpleSim) implemented in the Offline 
simulation framework. With this approach the lateral distribution of the air 
shower is assumed to follow a NKG-like functional form [29,30]. A model gen- 
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Fig. 5. Comparison between CORSIKA/Geant4 simulations and the fast 
CONEX/SdSimpleSim approach, (a): distribution of the time at which the SD 
station is triggered, (b): difference between the simulated and the reconstructed 
energies using the hybrid technique. The figures refer to events at logio(E/eV) = 
18.5. 

erating realistic signal timing for the closest station to the shower axis has 
been derived from a full Monte Carlo using AIRES [31] simulations. The Sd- 
SimpleSim code also includes the simulation of noise triggered stations, which 
could spoil the reconstruction of the event. The noise rate of the surface de- 
tector is self-adjusting to yield 20 Hz per station. As a cross-check, the number 
of noise triggered stations has been derived from data and the obtained dis- 
tributions have been parameterized. 

Dedicated CORSIKA/Geant4 simulations have been carried out to validate 
the performance of this fast approach against the full Monte Carlo method. 
Figure 5 shows the distribution of the station trigger times and the differ- 
ence between simulated and reconstructed energy as obtained with the two 
simulation modes. The consistency between these results provides a robust val- 
idation of the fast approach and makes it possible to produce of huge number 
of simulated events. 



5.3 Time Dependent Detector Simulation 

The Monte Carlo simulation for the calculation of the hybrid exposure has 
been based on the fast simulation approach described above. In fact for cov- 
ering all the energy ranges and the phase space of the detector configurations 
with enough statistical power, the number of simulated events is required to 
be largely oversampled with respect to the available raw data. The simulation 
has been designed to reproduce the actual sequencing of the detector status 
with a resolution of 10 min which corresponds to the time bin slot used for the 
on-time calculation. First a time is randomly chosen within the sidereal time 
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interval we want to simulate. Then all the relevant status information about 
each detector is retrieved from the on-time calculations. Based on the on-time 
fraction during the simulated time bin, only a sub-sample of the events is sent 
to the detector simulation. 

The CONEX showers used for this purpose have been generated from 10 17 
up to 10 21 eV. QGSJET-II [20,21] and Sibyll [32] have been used as high en- 
ergy interaction models. Proton and iron particles are taken as cosmic ray 
primaries. 

To account for the growth of the array with time and problems during the SD 
data-taking, only the active SD stations are considered during simulation. 

For the FD time dependent simulations the values of variance, baseline and 
trigger threshold averaged over lOmin are considered. The available FD ab- 
solute calibration data are used to adjust the simulated electronic gains on a 
pixel by pixel basis. This scales the shower signal with respect to the FADC 
trace noise and therefore influences the signal-to-noise ratio. In addition, in- 
correct cabling in some FD cameras is simulated for the instances discovered 
in the real detector. Data from the atmospheric monitoring system is used to 
set the hourly aerosol density profile as measured by the CLF [18] and the 
monthly mean molecular atmosphere as provided by balloon flights [33]. 



6 Event selection and validation of Monte Carlo simulation 

An unbiased measurement of the cosmic ray flux requires an exposure as free 
as possible from systematics. To this aim only high quality hybrid events are 
used. 



6.1 Quality cuts 

In this analysis high quality hybrid events have been selected using the fol- 
lowing criteria: 

• since we use the Gaisser-Hillas function [12] to evaluate the total calorimet- 
ric energy, a successful fit of the longitudinal profile with this function is 
required. Moreover, the y 2 per degree of freedom of the fitted profile should 
be less than 2.5; 

• the energy and shower maximum can only be reliably measured if X max is in 
the field of view (FOV) of the telescopes (covering 1.5° to 30° in elevation). 
Events for which only the rising or falling edge of the profile is detected are 
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Fig. 6. Relative difference between proton and iron exposure with respect to a mean 
composition exposure, as obtained from the Time Dependent Detector Simulation. 
In the left panel the dependence on the primary mass is clearly visible. In the right 
panel this difference is strongly reduced by the fiducial volume cut. 



not used, i.e. it is required that the depth of shower maximum be within 
the minimum and maximum observed depths; 

• to avoid potential systematic uncertainties related to the calculation of the 
Cherenkov light contribution, events with a relative amount of reconstructed 
Cherenkov light exceeding 50% of the total received light are not used in 
this analysis; 

• a good energy resolution is assured by accepting only events for which the 
total uncertainty of the reconstructed energy (including the propagated sta- 
tistical uncertainties of the detected photons, the geometry and the atmo- 
sphere) is smaller than 20%. 



Furthermore it is required that: 



• the aerosol content of the atmosphere is measured [18,34] for the time period 
of the event to allow a precise calculation of the transmission and scattering 
of photons in the atmosphere; 

• the cloud coverage according to Lidar measurements [34] is lower than 25% 
at the time of the event, since clouds could obscure part of the longitudinal 
profile and lead to an event selection inefficiency not accounted for in the 
aperture simulation. 



6.2 Fiducial cuts 



In addition to the above mentioned quality criteria, fiducial cuts have been 
applied to assure a robust calculation of the exposure, independent of the 
trigger threshold, mass composition and energy scale uncertainty. 
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To ensure that the probability of a trigger from at least one surface detector 
station is unity regardless of the primary particle, it is required that 



• the energy of the shower is larger than 10 18 eV; 

• the zenith angle of the shower is less than 60°; 

• the position of the station used for the hybrid reconstruction is within 
1500 m of the shower axis. 

The limited field of view of the fluorescence detector and the requirement of 
observing the shower maximum may both introduce a different selection effi- 
ciency for different primary masses. For instance, protons develop deeper into 
the atmosphere and have a deeper shower maximum than heavy primaries, on 
average. For vertical events the fraction of events with their maxima falling 
below the observation level is thus larger for proton primaries and correspond- 
ingly the selection efficiency is smaller. The mass dependence of the exposure 
for showers selected only by quality cuts is clearly visible in figure 6a. At low 
energies, where events are only detected close to the detector, iron primaries 
have a smaller exposure because of their shallower X max that tends to be more 
often above the upper limit of the FD field of view than it is for protons. At 
high energies the majority of the showers are far away from the telescopes and 
the bias is dominated by the lower field of view boundary that disfavors the 
selection of primary protons. 

In order to achieve an almost equal detection probability for all possible pri- 
maries, the following fiducial field of view cut has been designed: 

X up [g/cm 2 ] > 900 + 6 • (e - 18) (8) 

9 f 550 - 61 • (e - 19.06) 2 for e < 19.06 
X low [g/cm 2 ] < (9) 

550 for e > 19.06 



where e = \og 10 (E/ eV), and X up and Xi ow are the upper and lower boundaries 
of the telescope field of view which depend on the shower geometry. The 
application of this cut reduces the primary mass dependence above 10 18 eV. 
This is shown in figure 6b. 

A further cut was introduced in order to remove the FD trigger threshold 
effects induced by the energy scale uncertainties. The fluorescence detector 
trigger is in fact fully efficient for short distances between the shower and the 
detector. At larger distances the trigger probability decreases. A possible sys- 
tematic shift in the assignment of the shower primary energy, due to the 22% 
energy scale uncertainty (see tab. 1), may alter the derived trigger threshold 
and the exposure. To quantify this effect, the energy assignment of the simu- 
lated events can be shifted up and down by 22% before applying the selection 
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Fig. 7. The systematic uncertainty of the absolute energy scale of about 22 % would 
cause significant systematic uncertainties in the fiducial volume and the exposure 
(left panel). Dedicated event selection criteria are used to remove this dependency 
(right panel). The arrows show the cut values for the different energies based on eq. 
(10). 

criteria. The effect is shown in figure 7a. 

The possible dependence of the trigger threshold on a systematic shift in the 
energy assignment has been removed by dedicated selection criteria obtained 
from Monte Carlo studies. The available detection volume is limited by a set 
of fiducial volume cuts which require the shower core to lie within a distance 
-D max from the fluorescence detectors: 



L> max [km] < 



24 + 12(e- 19) 



for e < 18.5 



24 + 12(e-19) + 6(e-18.5) for e > 18.5 



(10) 



As is clearly shown in figure 7b this cut limits the available detection volume 
to a region in which the fluorescence trigger is saturated even if the energy 
scale is changed within the known systematic uncertainties (±22%). The expo- 
sure calculation thus becomes independent of the trigger threshold and of the 
absolute energy scale within current estimations of its systematic uncertainty. 



6.3 Cross checks 



All the above criteria have been applied to both data and MC events. The 
reliability of the quality criteria are checked by comparing the cut parameter 
distributions of data and Monte Carlo. Examples of these comparisons, shown 
in figure 8, indicate a very good agreement between data and Monte Carlo. 
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Fig. 8. Examples showing the agreement between simulation and data. Proton and 
iron primaries are shown separately for the simulated data. For each figure all quality 
and fiducial cuts are applied except the one related to the variable shown. The arrow 
denotes the selection cut on this variable. 



The exposure calculation depends somewhat on the hadronic interaction model 
used in the Monte Carlo simulation. Different hadronic interaction models pre- 
dict different fractions of shower energy converted into visible light [13] pro- 
ducing different energy assignments and X max predictions. These differences 
might affect the selection efficiency and lead to a model dependence in the ex- 
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Fig. 9. Relative difference between Sibyll and QGSJetll exposure with respect to 
the average of the two simulations. All quality and fiducial cuts are applied. 

posure. Two models, QGSJetII-03 and Sibyll 2.1 have been used as input for 
the Time Dependent Detector Simulation and the selection efficiencies have 
been compared. As is shown in figure 9 the effect is lower than 2% averaged 
over the whole energy range. For this reason the exposure has been calculated 
averaging the Monte Carlo samples simulated with the different interaction 
models. 



As mentioned above, the exposure is calculated as a function of reconstructed 
energy to correct for distortions of the spectrum introduced by the reconstruc- 
tion of the shower energy. It is well known that this correction depends on the 
initial assumptions of the true distribution (see eg. [35]), i.e. on the energy 
distribution of the generated events N(E) in Eq. (3). The exposure has been 
calculated using different distributions for N(E): power laws with three dif- 
ferent spectral indexes (7 = 2, 3, 3.5), a broken power law and the combined 
spectrum measurement from the Pierre Auger Observatory [7]. The ratio of 
the resulting exposure S(E rec ) to the undistorted one, £(E gen ), has been calcu- 
lated for the different cases. From this analysis, the choice of the input spectra 
used in the Monte Carlo results in a systematic uncertainty lower than 2%. 

The availability in the Pierre Auger Observatory of two independent detection 
techniques allows an overall validation of the Monte Carlo simulation chain 
with data. As shown in [4] the surface detector trigger is 100% efficient above 
a few EeV. Since SD data are unaffected by the distance to the FD-site, 
light attenuation or clouds, the FD trigger and selection efficiency can thus 
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Fig. 10. Conditional probability P(FD|SD). Comparison between measured and 
simulated efficiencies. The simulation has been performed using both proton and 
iron primaries. 

be measured directly from the data. A set of high quality SD showers have 
been selected during the time periods with at least one FD-site taking data. 
Given this set of Nsb showers, we count the number of events that had at 
least one triggered telescope, iV(FD tr i g ), and fulfilled all the selection criteria 
previously described, iV(FD se i). The FD trigger and selection efficiencies can 
then be estimated from: 



and 



£trig = P(FD trig |SD) 



e scl = P(FD scl |SD) 



iV(FD 



trig) 



SD 



iV(FD 



selJ 



iVc 



SD 



11^ 



(12) 



For each data shower, 20 simulated CONEX showers are generated with the 
given SD energy for proton and iron primaries. These showers are then pro- 
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Fig. 11. Hybrid exposure between November 2005 and May 2008 for proton and 
iron primary particles. 

cessed through the Time Dependent Detector Simulation with the same arrival 
time, direction and core position as measured by the SD, yielding the expected 
efficiencies: 



^ MC = PCFD.Igen) = *2™0 , (13) 

* *gen 

where i stands for either the trigger or selection criteria and iVg en denotes the 
number of generated events. 

The Monte Carlo prediction is compared with the measurements in figure 10. 
It can be seen that the shape of the two curves agree for both proton and iron 
simulations. However we note a normalization factor between simulation and 
data of 0.92 ± 0.02 assuming a mixed composition of 50% proton and 50% 
iron nuclei. This could be related to an uncertainty in the on-time, or caused 
by the poorer energy resolution of the SD. 



7 Exposure 



The hybrid exposure is shown in figure 11 for both proton (full circles) and 
iron (open squares) primaries. It is calculated for the data period between 
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November 2005 and May 2008, and is that used for the hybrid energy spectrum 
measurement published in [7]. The analysis of the Central Laser Facility shots 
described in section 4 has revealed a systematic shift in the on-time calculation. 
To take account of this effect the exposure has been reduced by 4%. Moreover 
the end-to-end comparison in section 6.3 has shown that the ratio of the true 
event rate to that expected from Monte Carlo is 0.92 ± 0.02. The systematic 
uncertainty of this comparison has been estimated to be ± 5%. Consequently 
the exposure has been reduced by half of the corresponding correction (~ 4%) 
to cover the full range of expectations. These two corrections are included in 
the exposure shown in figure 11. 

A mixed composition of 50 % proton and 50 % iron nuclei has been assumed in 
the exposure calculation [7]. Numerical values of the hybrid exposure can be 
found in [36]. The remaining composition dependence has been included in the 
systematic uncertainty. This was found to be about 8% at 10 18 eV decreasing 
down to 1% above 10 19 eV (figure 6b). The dependence of the exposure on 
the hadronic interaction model has been studied in section 6.3. The effect 
is smaller than 2% over the entire energy range used for the calculation of 
the exposure. The dependence of the exposure on the different input spectra 
used in the Monte Carlo simulation has also been investigated and found to 
be smaller than 2%. The overall systematic uncertainty on our knowledge of 
the hybrid exposure has been obtained by summing all these contributions in 
quadrature. It ranges from about 10% at 10 18 eV to 6% above 10 19 eV. 

In figure 12, the growth of the hybrid exposure as a function of time is shown 
for three different energies. The increase with time shown at each energy comes 
as a result of the concurrence of different effects, i.e the accumulation of data 
taking with time and the growth of the SD array. One can also observe faster 
changes corresponding to the longer FD data-taking periods in the austral 
winter. The effect due to the growth of the SD array is more marked at higher 
energies where a larger hybrid detection volume is accessible with the new SD 
stations. 



8 Summary 

A method for the calculation of the hybrid exposure of the Pierre Auger Obser- 
vatory has been developed. The method is mainly based on the Time Depen- 
dent Monte Carlo simulation technique. This technique allows the simulation 
of a sample of events that reproduces in detail the exact configuration of the 
experimental data taking, including both instrumental and atmospheric con- 
ditions. 

With this aim the on-time of the hybrid detector has been calculated in a 
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Fig. 12. The growth of the hybrid exposure as a function of time starting from April 
2006 up to May 2008 for three different energies. 

very accurate way taking into account environmental and instrumental effects 
occurring at different levels of the DAQ process, with respect to both the FD 
and the SD. The on-time information has then been used as input for the 
Time Dependent Monte Carlo simulation. 



In order to follow the fast changes of the detector configuration and to get 
acceptable statistical and systematic uncertainties, a simulation with very 
high statistics is crucial. A fast simulation has been used, using CONEX 
shower profiles for the FD simulation and an efficient simulation of the SD 
response. This fast approach has been validated using a dedicated CORSIKA 
plus Geant4 simulation. No significant difference has been found between the 
two approaches. 



To obtain an unbiased measurement of the cosmic ray flux the exposure es- 
timate must be as free as possible of systematics. For this reason only very 
high quality events have been used [7]. To satisfy this aim a set of quality 
criteria has been developed and discussed in the paper. The effect of the qual- 
ity criteria has been cross-checked by comparing the data and Monte Carlo 
distributions; a very good agreement has been found. 

The systematic uncertainties arising from the unknown details of mass com- 
position, hadronic interaction physics and the true energy spectrum have been 
calculated, and residual systematics from the on-time calculation have been 
estimated. The overall systematic uncertainty on the calculation of the hybrid 
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exposure has been found to be lower than 10% (6%) at 10 18 eV (above 10 19 eV). 
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